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Abstract 

The classical theory of radiating point-charges is revisited: the retarded 
potentials, fields, and currents are defined as nonlinear generalized func- 
tions. All calculations are made in a Colombeau algebra, and the spinor 
representations provided by the biquaternion formulation of classical elec- 
trodynamics are used to make all four-dimensional integrations exactly and 
in closed-form. The total rate of energy-momentum radiated by an arbi- 
trarily moving relativistic point-charge under the effect of its own field is 
shown to be rigorously equal to minus the self-interaction force due to that 
field. This solves, without changing anything in Maxwell's theory, numerous 
long-standing problems going back to more than a century. As an imme- 
diate application an unambiguous derivation of the Lorentz-Dirac equation 
of motion is given, and the origin of the problem with the Schott term is 
explained: it was due to the fact that the correct self-energy of a point charge 
is not the Coulomb self-energy, but an integral over a delta-squared function 
which yields a finite contribution to the Schott term that is either absent or 
incorrect in the customary formulations. 

03.50.De Classical electromagnetism, Maxwell equations 
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1 Introduction 



This is the long version of the paper The self-interaction force on an arbitrarily 
moving point-charge and its energy -momentum radiation rate: A mathematically 
rigorous derivation of the Lorentz-Dirac equation of motion, Ref. 0], in which 
the results and the main steps in the calculations were given in tensor form. In 
the present paper we give the details of these calculations, which because of their 
complexity are made using Hamilton's biquaternion formalism^ 

Gauss's theorem in its various forms is a major tool in field theory because it 
enables to transform integrals over closed surfaces into integrals over the volume 
contained within them. This enables, for example, to express Maxwell's equations 
in integral rather than differential form, and thus to establish precise relations 
between the total flux of fundamental field quantities through closed surfaces and 
the volume integrals of quantities that are differentially related to them. 

As long as all field quantities are smooth enough, and provided all surfaces are 
sufficiently regular and simple, there are no problems: Gauss's theorem and its 
numerous corollaries can be applied without restriction and with no great technical 
difficulties. But this is not the case when there are singularities in the fields. 

Consider for example the flux of energy-momentum radiated by a point-charge 
under the effect of its own field through a 3-surface enclosing its world-line, and 
the related 4- volume integral of the self-interaction force on this point-charge due 
to that field: according to Gauss's theorem they should be equal — provided all 
problems related to singularities are properly taken care of since the fields are 
infinite on the world-line. Indeed, let T( ) by the symmetric energy-momentum 
tensor, and 3 S the 3-surface enclosing a 4-volume 4 f2 surrounding the world-line 
between the proper-times t\ and r 2 . The radiated energy-momentum flux is then 

P(£):=/// T(d 3 S)= r dir III rf 3 fiT(V), (1.1) 

where Gauss's theorem was used to write the right-hand side. Using the definition 
of T( ) in terms of the electromagnetic field bivector F = E + iB, i.e., 

T() = ^-F+()F, (1.2) 

and Maxwell's equations in the form 

VF = -4vrJ, which implies F + V = +4vrJ, (1.3) 



'To get the Minkowski metric in that formalism all time-like scalars are pure imaginary. Thus 
the proper time is it, the proper-time derivative is ( ) := dQ/dir, and the proper-time integral is 
J ( )dir. For an introduction to Hamilton's quaternions see 0. 
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where J = J + is the charge-current density, the right-hand side of Eq. (11.11) 
becomes 

~ / 2 dir JJL d3n j ^ := / 2 dir = (L4) 

Therefore, comparing with the left-hand side of Eq. (11.11) . Gauss's theorem implies 

P(E) = P(Q) = -Q{Q), (1.5) 

where 

Q(Sl):= JJJ^ d 3 n^(F + J -JF). (1.6) 

Equation Eq. (11.51) states that the rate of energy-momentum flux through the 
3 -surface 3 S at the time r, calculated as the proper-time derivatives of either 
the 3-surface integral P(S) or of the 4-volume integral P(f2), is equal to minus 
the integral Q(tt) of the self-interaction force density \{F + J — JF) over the 
3-volume 3 f2 enclosed by the 3-surface at that time. 

However, when -P(S) and P(Cl) are evaluated using for F the customary 
Lienard-Wiechert field strengths and for J the customary current density of an 
arbitrarily moving relativistic point-charge, what is not without difficulties and 
ambiguities since both F and J are unbounded on the world-line, one finds that 
they are not equal. 

For example, following Lorentz [3], the self-force is usually calculated non- 
relativistically and then generalized to the relativistic case. This gives 

P(Q) = lim ~Z - \e 2 iA 2 Z - \e 2 %Z\ (1.7) 

where the term containing the square of the acceleration A 2 = \A\ 2 = ZZ, 
known as the 'Larmor term,' is the invariant radiation rate, and the term in the 
third proper-time derivative of the position, Lorentz's frictional force, is known as 
the 'Schott term.' 

On the other hand, the energy-momentum flux rate can be found by either 
generalizing a non-relativistic calculation, or, following Dirac JH, by a direct 
covariant 4-dimensional calculation. This gives 

P(E) = lim^Z-^e 2 M 2 Z, (1.8) 

in which the Schott term is missing and where the diverging quantity e 2 /2£ is 
called the 'Coulomb self-energy' or 'electrostatic self-mass' of the point-charge, 
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which differs from the 'electromagnetic self-mass' 2e 2 /3£ which appears in (11.71) 
by a factor 4/3. 

Apart from this notorious '4/3 problem,' the main difference between -P(S) 
and P(fi) is that is a Minkowski force, that is a genuine relativistic force 
satisfying the orthogonality condition 

PoZ = 0, (1.9) 

on account of the kinematical identities ZZ — 1, Z o Z — 0, and Z o Z — —A 2 , 
whereas -P(S) is not such a force because there is no Schott term. In fact, this 
missing term is a major problem because, as noted by Dirac, the derivation of 
(|1.8I) is quite simple, and there is no way to get the Schott term in this calculation 
without making a supplementary assumption. 

These problems, as well as other inconsistencies, have led to hundreds of 
papers and endless controversies which we are not going to review and discuss 
here. Indeed, while theses difficulties are known, there is a large consensus that 
the finite part of the force (11.71) is the correct relativistic form of the self-force 
on an arbitrarily moving point-charge, i.e., the so-called 'Abraham' or 'von Laue' 
self-interaction force |ie 2 (Z + A 2 Z). Similarly, there is a large consensus that 
the divergent terms in (11.81) and (11.71) should be discarded or 'renormalized', and 
that the equation motion of a point-charge in an external electromagnetic field F ext 
is 

mZ = ^ie 2 (z + A 2 Z) + e^{F+ t Z - ZF ext ), (1.10) 
i.e., the Abraham-Lorentz-Dirac equation of motion. 

In this paper we show that the fundamental reason for discrepancies such as 
between (11.81) and (11.71) is mathematical, and due to two shortcomings of the 
customary calculations: 

• incomplete characterization of the singularities of the fields and current at 
the position of the point charge — everything should be defined not just as 
distributions but as nonlinear generalized functions ; 

• use of formalisms such as elementary calculus or Schwartz distribution 
theory in which the product of distributions is not defined — this is a major 
defect since the energy-momentum tensor (|1.2I) is nonlinear in the fields so 
that all calculations should be made in a Colombeau algebra imi8ll9l[i~0l. 

Once all singularities are properly defined, what basically consists of redefining 
the Lienard-Wiechert potential as a Colombeau nonlinear generalized function 
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[fTTl[T2l . the calculations of the self-energy-momentum radiation rate and of the 
self-interaction force are straightforward (but technically complicated) and the 
equality (11.51) is obtained for a strictly point-like electric charge in fully arbitrary 
relativistic motion without making any modification to Maxwell's theory. 

The importance of this result lies is the fact that it leads to the solution of 
a number of difficulties which have often been interpreted as suggesting that 
Maxwell's electrodynamics of point-like electrons is not fully consistent, and 
required to be modified at the fundamental level. To illustrate that this is not 
necessary, it will be shown how the Lorentz-Dirac equation of motion derives 
from our rigorously calculated self-interaction force. 

Because the full calculation leading to the proof of the equality (11.51) for a 
relativistic point-charge is very long and complicated, the details were not given 
in the paper 0, where only the general physical concepts, the principle of the 
mathematical methods used, and the main steps in the calculation were given. 



2 Definitions of the fields and of the geometry 

Let Z be the 4-position of an arbitrarily moving relativistic point-charge, and X 
a point of observation. At the point X the 4-potential A, field F, and 4-current 
density J are functions of the null interval R between X and Y, i.e., 

R := X - Z, such that RR = 0, (2.1) 

as well as of the 4-velocity Z, 4- acceleration Z, and 4-biacceleration Z of the 
charge, to which three invariants are associated: Z o /?, Z o R, and Z ° R. The 
first one is called the retarded distance, 

i£:= Zo R, (2.2) 

which enables to introduce a 'unit' null 4-vector K defined as 

K(6,4>):=R/Z, (2.3) 

and the so-called acceleration and biacceleration invariants defined as 

k:= Z o~K, and %x ■= Z ° ~K. (2.4) 

The three spatial variables {£, 9, <f\ and the proper time it define a causal 
coordinate system that is most appropriate for expressing the fields of a point- 
charge moving according to the world-line prescribed by Z(r). As explained in 
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Refs. IfTTTl and lfT2ll . the retarded potential of such a point-charge must be written 
as a nonlinear generalized function, i.e., 



(2.5) 



where T(£), which is absent in the customary Lienard-Wiechert formulation, is 
the generalized function defined as 



undefined r < 0, 



T(r) := { r = 0, and -^T(r) = 5(r), (2.6) 



r > 0, 



which explicitly specifies how to consistently differentiate at £ = 0. Instead of the 
customary Lienard-Wiechert field, the field strength is then lfT2l 



-e((l-^-^)T( 
+ e((l- 



(2.7) 



which apart from the presence of the T-function factor, has an additional 5-like 
contribution. As a consequence of these T and 5 factors, the current density 
deriving from this field, i.e., 
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r=r r 



(2.8) 



is significantly more complicated than the customary current density, but turns 
out to be locally conserved, whereas the customary one is not [|TT1 [T2| . In this 
equation, as in Eqs. (12.61) and (12.71) . the condition r = r r implies that all quantities 
are evaluated at the retarded proper time r r . In the following, for simplicity, this 
condition will be specified explicitly only for the main equations. 

Having defined the field F and the current density J to be used to calculate 
the integrands of Eq. (|l.ll) . (11.41) and (11.61) . it remains to define the geometry of a 
3-surface enclosing the world-line in order to calculates the integrals themselves. 

In view of this we begin by ignoring the fact that F is singular on the world- 
line. We are then free to chose any closed 3-surface intersecting the world-line 
at two fixed proper times T\ and r 2 . This is because by Gauss's theorem the 
field-integrals over the 4-volume enclosed between any two such 3-surfaces will 
give zero since VF = 0, and thus T(V) = 0, in that 4-volume. To be consistent 
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with the already defined causal coordinate system {ir, £, 6, <p} we therefore take a 
3-surface consisting of one 3-tube of constant retarded distance £ = Cst between 
the times n and r 2 , and, to close that 3-tube at these times, two 3-spheres r = Cst 
at the times n and r 2 between the world-line and the radius of the 3-tube. We 
therefore do not take a light-cone (defined as r r := r — £ = Cst) to close the tube 
at both ends as is often done to simplify calculations when only surface integrals 
are calculated: such a choice would make the corresponding 4-volume integrals 
significantly more complicated. Besides, as the formulation of our fields involves 
T and 5 functions whose argument is £, taking surfaces and volumes different 
from those defined by constant values of the causal coordinates £ and r would 
make the evaluation of the fields singularities at £ = even more complicated 
than they already turn out to be. Finally, with our choice, the 4-volume element is 
the simplest possible, namely 

d 4 n = dir d(f) d6 sin(0) d£ £ 2 , (2.9) 

so that the 4-volume integrals are simply proper-time integrals times ordinary 

3- volume integrals with the radial variable equal to £. 

Now, to deal with the singularities on the world-line, we replace the just defined 

4- volume, that is such that £ varies between £ = and some maximum £ = £ 2 , 
by a 'hollow' 4-volume in which £ varies between £ = £i ^ and the maximum 
£ = £ 2 . This 4-volume will then be comprised between two 3-tubes of radii £l and 
£ 2 , where £i will eventually tend to at the end of the calculation. Thus, we will 
calculate the 4-volume integral 

T2 rrr / \ r T2 r& 

2 



dir JJJ d 6 Q (^...J = J dir J £ JJ duJ [■■■), ( 2 -!0) 

in which (for practical reasons) the integration over the ordinary solid angle duj = 
d(j> d9 sin(6 ) ) will be done first in general. 

In this geometry, the normal to the 3-tube £ = Cst is§ 



iV T = V£ = (1 - k£)K - iZ 
and the integrals over the internal and external 3-tubes are 

' d?Y? (...) =J T2 dir jj due(-)N T 



(2.11) 



(2.12) 



2 In the causal coordinate system {ir, £, 8, </)} the operator V is the causal 4-differentiation 
operator defined in [12, Sec. 2.3]. 
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where ^ is the internal radius of the tube that is sent to zero at the end of the 
calculation. 



Similarly, the normal to the 3-sphere ir = Cst is 

N s = Vir = -in^K + Z 



(2.13) 



and the combined integrals over the two hollow 3-spheres enclosing the hollow 
4-volume at both ends can be written 



(2.14) 



3 Mathematical methodology 

The proof of Eq. (11.51) involves calculating products of the generalized functions 
T, S, and 5' such as T 2 , ST, S'T, S'S, and S 2 , possibly multiplied by £ n with 
n G Z. If these generalized functions were defined as Schwartz distributions these 
products would not be defined and all sorts of inconsistencies would arise. But as 
we define them as nonlinear generalized functions we can use Colombeau's theory 
which enables to multiply, differentiate, and integrate them as if we were working 
with ordinary smooth functions f71l8ll9l [TOll . 

In Colombeau's theory the generalized functions T, 5, etc., are embedded in 
an algebra by convolution with a 'mollifier function' r/(x), i.e., an unspecified real 
smooth function whose minimal properties are given by 

Jdzr){z) = l, and J dz z n r](z) = 0, Vn = 1, q G N. (3.1) 

This function can be seen as describing details of the singularities at the infinites- 
imal level. These details are possibly unknown, but not essential for getting 
meaningful results which are generally obtained by integrations in which only 
global properties such as (|3.1I) are important. However, as will be the case in 
this paper, further constraints can be put on the mollifier, for instance to meet 
requirements imposed by the physical problem investigated. 

Of course, working in a Colombeau algebra requires some precautions: all 
'infinitesimal' terms aring as a result of a differentiation or any other operation must 
be kept until the end of the calculation — contrary to calculating with distributions 
where such terms may possibly be discarded. Similarly, distributional identities 
like 5'(x) = —6(x)/x, which was used to write the conserved Lienard-Wiechert 
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current density in the form (|2.8I) . have to be avoided to remain as general as 
possible. We will therefore use J in the form lfT2l Eq. (4.17)] 

where iV T is the normal to the 3-tube £ = Cst, i.e., Eq. (12.1 II) . 

To calculate the integrals (12.101) . (12.121) . and (|2.14l) . we introduce a notation 
that will be systematically used in the following. That is, we rewrite the field 
strength as a sum of two terms 

F(t,C) =F r (T,0T(0 + F 6 (T,0S(0, (3-3) 

where the upper indices ' T ' and l<5 ' refer to the T part (i.e., the customary Lienard- 
Wiechert field) and respectively to the 5 part of Eq. (12.71) . The energy-momentum 
tensor T() consists then of four terms, and consequently the surface integrals 
(12121) and (l2~T4l) comprise also four terms, i.e., P T,T (£), P T,5 (X), P 5,T (X), and 
P S ' S (E), in which the factors T 2 , TS, ST, and 5 2 appear. Similarly, the volume 
integral (|2.10l) involves then two terms P T ' 5 (fi) and P s ' s (Cl), in which the factors 
5T and 5'T, and respectively 55' and 5 2 , appear. 

As for the integrations themselves, we perform them according to a similar 
pattern, in which all £ integrations are left for the end because both 5 and T are 
functions of £ which are singular as £ — ► 0. Thus we begin by making the angular 
integrations, which together with the complicated algebra implied by the products 
of the many field components is the most laborious part of the calculation. This 
work is however greatly facilitated by the biquaternion techniques due to Paul 
Weiss explained in Sec.|4] which take maximum advantage of the cancellations 
between null 4- vectors, and which allow a maximal separation of variables and 
factorization of all quantities written in spinor form in terms of factors depending 
solely of kinematical or angular variables lfT3l . The resulting surface and volume 
integrals are given in Sees. [5] and [6] 

Next we transform the resulting surface integrals in such a way that they have a 
form directly comparable to the volume integrals, namely a double integration over 
£ and r. This enables, after dropping the r integration, to compare the remaining 
^-integrals, which should be equal according to Eq. (11.51) . The integrands of 
Eqs. (12.121) and (12.141) . which are initially integrated over r and respectively £ 
alone, have therefore to be transformed and combined in non-trivial ways. This is 
made possible by the fact that all kinematical quantities are implicit functions of 
the retarded time r r = r — £, so that their r and £ derivatives are related as 

— E(r r ) = -T^E(r r ), or E(r r ) = iE{r r )' ', (3.4) 
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where the dot and the prime denote differentiation with respect to ir and to £, 
respectively, and where E(r r ) is any expression function of the variable r r . 

A particularly interesting intermediate result in the transformation of the sur- 
face integrals is given as Eqs. (I7.9H7.111) : it shows that the Coulomb self-energy 
does not contribute to the total surface integral P(£). 

Finally, the volume integrals given in Sec. [6] are transformed in such a way 
that there is no 5' anymore. This is done using the identities (5 2 )' = 255' and 
(f 5T)' = 5T + £5'T + £5 2 , which enable to integrate (RT25b and (l6\26l) by parts 
to eliminate 5'. 

At this stage all expressions are comparable. It is then found that Eq. (11.11) 
holds exactly as it should according to Gauss's theorem, and that both P(X) and 
P(fi) can be written in the form (I7.3H7.6I) . Thus Eq. (11.51) holds as well and 
P(£) = P(f2) can be written in the form (I8.1H8.4I) . However, since we have 
not yet made the final £ integrations, these equalities are only a check that all 
algebraic calculations and angular integrations are correct, and that the integrands 
have been properly transformed using Eq. (13.41) and the rules of partial integration. 
The precautions taken to properly deal with the singularities, and the postulate 
that everything is calculated in a Colombeau algebra, will become essential when 
making these £ integrations and taking the limits £ 2 — ► oo and £x — > 0, what will 
be done in Sec. [8] 



4 Spinor representations 

So far we have used biquaternions as an alternative to the tensor notation. We now 
give the reason for using biquaternions in the present problem: it is to provide 
concise and explicit representations of all retarded quantities. That will greatly 
simplify the calculations and will lead to final results that are exact and in closed- 
form. 

The starting point is the spinor decomposition of the kinematical quantities 
(2j|. For instance, since biquaternions form an algebra, the 4- velocity has a unique 
decomposition as a product 

Z(t) = BB + , such that B = B + , and BB = 1. (4.1) 

This is the spinor representation (or decomposition) of the 4-velocity in the bi- 
quaternion algebra, in which the unit biquaternion £>(t) biunivocally defines (up 

3 Note that we write B + even though algebraically B + = B. This is to keep track of Lorentz 
covariance, which will be explicit in all our formulas. 



10 



to a factor ±1) the general Lorentz-boost to the frame in which the arbitrarily 
moving particle at the point Z(r) is instantaneously at rest@ Taking the imaginary 
proper-time derivative of both sides of Z = BB + we get 

Z{t) = BB + - BB + , (4.2) 

so that in the biquaternion algebra the 4- acceleration has the spinor representation 

Z( T ) = BAB + , where A = BB-B + B*, (4.3) 

where .4 is minus^l the invariant acceleration in the instantaneous rest- frame. 

To represent retarded quantities we need explicit covariant expressions for null 
4-vectors such as the causal interval (|2.1I) between the points X and Z. These 
were discovered in 1941 by Paul Weiss, who found that the spinor decomposition 
of the null interval X — Z is lfl"3l . 

X - Z = R = £B(i + v)B + . (4.4) 

Here z7 = P(9, 0) a unit vector so that i + V, and consequently R, are null 
biquaternions. Equation (|4.4t is therefore an explicit parametrization of R in 
terms of four variables: The invariants £ and r, and the two angles 9 and 
characterizing the unit vector V in the instantaneous rest frame of the particle. 
Thus, in that frame, i.e., when B — 1, the interval X — Z reduces to the quaternion 
£(i + V), which shows that in that frame the vector £z7 corresponds to the ordinary 
radius vector x — z, and that the distance £ appears to an observer at rest in that 
frame as the ordinary distance \x — z\. On the other hand, when B ^ 1, Eq. (14.41) 
provides a general parametrization of the null-interval R, that is, in geometrical 
language, of the 'light-' or null-cone originating from Z. 

In the case where Z, £, i>(9, 0), and B are evaluated at the retarded proper- 
time, Eq. (14.41) defines a causal coordinate system erected at the retarded 4-position 
Z(r r ), i.e., 

(4.5) 



x{T,z,e,<i>) = z + zB{i + v& 



Paul Weiss called it 'retarded coordinate system,' but we prefer the term causal 



coordinate system. From the form of equations (|4.4I) and (14.51) it is immediately 
seen that £ can be expressed in terms of R and Z, namely as the scalar product 



ii = ZoR 



(4.6) 



4 As explained in Ref. [2|, any unit biquaternion C, i.e., such that CC = 1, defines a proper 
Lorentz transformation. Then, a 4-vector is a biquaternion like Z which transforms as CZC + in a 
Lorentz transformation, whereas a spinor is a biquaternion like B which transforms as CB in such 
a transformation. 

5 This sign convention is used throughout the paper. 
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which confirms the usual definition (12.21) of the retarded distance £. Also, as 
K — R/£, and with the 4- acceleration expressed as (14.21) . it is easily found that 
the acceleration invariant (12.41) is 



k = V ■ A 



i 

r=T r 



(4.7) 



whereas there is no similar explicit formula for the biacceleration invariant \- 

For convenience, it is useful to write the null biquaternion i + V in terms of an 
idempotent a such that aa = a and aa = 0. Therefore, we write 

i + u = 2ia, where a(6, 4>) := ~(1 + Of), (4.8) 

so that R and the 'unit' null 4- vector K defined by (12.31) become 

R = 2i£BaB + 1 and K = 2iBaB + . (4.9) 

Inspecting the field, currents, and normals, i.e., Eqs. (12.71) . (12.81) . (12.111) . (12.131) . 



and (13.21) . we see that we have all we need to express all the quantities appearing 
in our problem in spinor form, namely by means of the spinor representations of 
Z, Z, and K, given by (|4.1I) . (|4.3I) . and (14.91) . For instance, after some simple 
algebra, the retarded field (12.71) can be put into the form 



so that 



^&*{v + 2ti*Aa)B + r{0 - |(1 - &.A)B*VB+5{Z), (4.10) 



F + = -^-B(u + 2^aAa)BT{0 + Ul - &-A)BVB5{£). (4.11) 

Similarly, the normals (12.1 II) and (12.131) can be written 

N T = B{p - 2£aAa)B + , (4.12) 
N s = B{1 - 2i£aAa)B + , (4.13) 

where A and B are functions of r r = r — £ only, and and a are functions of 
and only@ 

By various transformations Eqs. (14.101) to (14.131) could be put in other equiv- 
alent forms, for example by not replacing the invariant «(r r , 9, <p) by its explicit 



6 In the non-relativistic limit (i.e., v <C c) the normal to the 3-sphere becomes N^ R = 1 — 
2i£7fAo : 7^ 1. It is only in the Gallilean limit (i.e., c — > oo) that the 3-sphere becomes the 
'now-plane' whose normal Nq — 1 has all vector (i.e., space) components equal to zero. 
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expression (14.71) . However, the technical advantage of the formulations (14.91) to 
(|4.13l) is that they provide a maximal separation of variables and a factorization 
in terms of factors depending solely of kinematical or angular variables. This 
will enable to use identities such as BB = B*B + = 1 efficiently, and to integrate 
over the angular variables independently from the others. Moreover, since the null 
components of all quantities are expressed in terms of the idempotent a, there will 
be simplifications and cancellations when multiplying them due to the identities 
oo = a and aa = 0, as well as vo = —io and Va = ia, etc. 

In practice, when calculating the energy-momentum and self-interaction force 
integrals, one will reach a stage where, after having made all algebraic simpli- 
fications, it will be necessary to integrate over angles the resulting expressions 
containing products of several a, a, V, A, and v ■ A factors in various positions. 
For example, 

1 ff . r_ 1, . I ff. 1 



,. dtu a A a = -A , whereas — // dtu a A a = -A. (4.14) 
47r JJ 3 4vr JJ 6 

The corresponding integrations are in general elementary, but they become easily 
laborious when the number of factors is larger than three. The result of a number 
of such integrations is given in Table [T] 

A special problem arises when integrating expressions containing the biaccel- 



eration invariant x = —iZ ° K, which appears in the current densities (12.81) and 



(|3.2I) . because there is no simple spinor representation for Z comparable to (14.31) 
for Z. Nevertheless, referring to (|4.4I) or (|4.9I) . one can write 

K = iZ + S, where S := BvB + , (4.15) 

and therefore, because of the kinematical identity Z ° Z = —A 2 , 

X = -iZ oK = -A 2 -iZ oS. (4.16) 



5 Surface integral: the energy-momentum flow 

In this section we calculate the four contributions P T T , P T > 5 , P 5 T , and P s,s to 
the surface integral -P(X), i.e., Eq. (11.11) or the proper-time integral of the left-hand 
side of (11.51) . The first one will be calculated in full detail to illustrate the method. 
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5.1 Surface integral P 



T,T 



We start from the definitions (ll.H - ll.2l) and calculate first the integral over the 
3-tube S T , i.e., 



P 



T (d 3 S T ) = — /// F r+ d 3 E T F r 
Tj T 87r JJJtt 



in 



dir // dcu £ 2 F r+ N T F r 



6 



(5.1) 



where we used (|2.12l) for integrating. With the T part of the field taken from 
ti4.10J - l4~TTl) . and the normal to the 3-tube given by (14.121) . this becomes 



P T ' T (E T ) = / dir II duj( 2 —B(-v-2£aAa)BT(i) 



in 



x B{y- 2£aAa)B + — B*(u + 2£aAa)B + T(£) 



Using the kinematical identities BB = B*B + = 1 and rearranging we get 

x B(-v- 2£aAo)(V- 2£oAo){v + 2£aAa)B + 
which after multiplying the first two parentheses becomes 



x B{1 + 2i£oA){v + 2£aAo)B^ 



6 



and then 



P 



87T 



^ ) = - dir dco -T%) 



x B [P+ 2&A W + 2i£aAP + Ai^WA aAa)B 



&2 



The angular integrals are then made using the formulas in Tabled] to get 



Ti 



x B( + 2^A + 2^A + &e{- l -A 2 ) )& 



(5.2) 



(5.3) 



(5.4) 



(5.5) 



(5.6) 
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which using (|4.1I) and (14.31) to revert from spinor into vector notation finally gives 
P T ' T (S T ) = e 2 £ rfzr T 2 (0(^ - \iA 2 z). (5.7) 



We now repeat the same procedure for the integrals over the two 3-spheres T, s 
closing the 3-tube at both ends, i.e., 



p T - T (s 5 ) = yyy r(<f£ 
i S r6 



57T 



F T+ rf 3 S T F T 



— J <li II (U- £ 2 F r+ N s F r 



(5.8) 



where we used (|2.14l) for integrating. With the T part of the field taken from 
ti4.10J - l4~TTl) . and the normal to the 3-spheres given by (|4.13l) . this becomes 



P T - T (S S ) 



8tt 



6 



6 



x B(l - 2i^aAa)B + —B*(u + 2^aAa)B + T^) 



(5.9) 



Using the kinematical identities BB = B*B + = 1 and rearranging we get 



x £(-i7- 2fov£a)(l - 2i£oAa){v + 2^a^a)i3 + 



(5.10) 



which after multiplying the first two parentheses becomes 



i y & /•/■ , e 2 



•57T 



6 



<K II doj -T 2 (0 



x B{-V - 2£aA) {V + 2£ala)B^ 



T-2 



(5.11) 



and then 



in 



6 



£ 2 



x ,6(1 + 2«£<j.Aa - 2^aAv - ^ 2 aAaAa)B 



(5.12) 
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The angular integrals are then made using the formulas in Tabled] to get 



• B[ 1 + 2i^A + 2z£jU - 4f{~A 2 ))B" 



(5.13) 

which using (|4.1I) and (14.31) to revert from spinor into vector notation finally gives 

^2 / 1 . 1 .. 9 . \ T 2 



P T > T f£ T ) 



Si 



I • I •• 2 

— Z + — iZ + -A 2 Z 
2£ 2 2£ 3 



(5.14) 



Therefore, adding (15.71) and (15.141) . the p T > T energy-momentum flux through 
the closed hypersurface is 

6 



P T ' T (£) = e 



2 i T2 dzr (±z-hA 2 z)r 2 



Tl 

6 



fi 



ei 



2f 



Z + — iZ + 
2£ 



U 2 Z)T' 



T2 



(5.15) 
(5.16) 



5.2 Surface integral P s ,T and P 



)T,<5 



Before calculating the mixed contributions P 5 ' r and P T ' S we remark that energy- 
momentum is a bireal quantity, i.e., P = P + . Thus if either P s ' r or P r ' 5 is 
calculated and found to be bireal, P 5 ' r and P r ' s will necessarily be equal to each 
other because 



P 5 ' T (E) 



87T 
1 

8^ 



F 5+ d 3 J:F r 



P T+ d 3 £P d = P' 1,d (E) 



T,<5/ 



Thus we begin by calculating the integral P s ' r over the 3-tube £ T , i.e., 



P 5 ' T fE T ) 



•57T 



rfir jj dtu £ 2 F d+ N T F 



2 AT T I? T 



(5.17) 



(5.18) 



With the 5 and T part of the field taken from (14.101 - 14311) . and the normal to the 
3-tube given by (14.121) . this becomes after some simplifications 

„2 



P 5 > T (E T ) = - / dir II du ^<5(0T(0 



in 



x B(l - iv ■ A)V(V - 2£aAa)(P + 2£aAa)B^ 



(5.19) 
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which after some algebra gives 

I r T 2 rr „i 



p<5,T (s T^ 



in 



xB[-P- 2i£PAa + £{v ■ A)P + 2ii 2 {y ■ A)VAo ) B 



£ 2 



6 



. (5.20) 



The angular integrals are then made using the formulas in Tabled] to get 



P 



1 f T2 e 2 
= 2j dlT e mT{0 

x6(0-^I + ^I + ^ 2 ))8 H 



6 



(5.21) 



which using (14.11) and (14.31) to revert from spinor into vector notation finally gives 

£2 



P 



<5,T /vVT\ 



z 2 £diT5(Z)r(Z)(-~iA 2 Z 



£1 



(5.22) 



We now repeat the same procedure for the integrals over the two 3-spheres E s 
closing the 3-tube at both ends, i.e., 



P ^(E S ) = — I dUI dueF d+ N»F 



1 



57T 



£2 



£1 



-2 t?8+ t\tS zpT 



(5.23) 



With the 5 and T part of the field taken from d4.10M4TTTT) . and the normal to the 
3-sphere given by (14.131) . this becomes after some simplifications 



x B(l - iv- A)V{\ - 2i£aAa)(P + 2&Aa)B^ 



T-2 



(5.24) 



which after some algebra gives 

r& rr . 2 



pS ' r ^ s ) = ^-J^dZ J! du e - 2 5(OT(0 

x B(-l + 2£PAa + £(P-A) -2S, 2 (v ■ A)vAa)B + 



(5.25) 



n 



The angular integrals are then made using the formulas in Table Q] to get 



1 

4 



P 5 ' T (S 5 ) = - ^ d£ ^5(0T(0 



xB(-i-ze^+o-e 2 (-^ 2 ))^ 



(5.26) 
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which using (|4.1I) and (14.31) to revert from spinor into vector notation finally gives 



P 



6 



6 



6' 



(5.27) 



As was explained at the beginning of that subsection, if P S ' T is bireal, what is 
the case of (15.221) and (15.271) . it can be taken as equal to P r ' 5 . Therefore, doubling 
the sum of (15.221) and (15.271) . the P r ' 5 + P 5 ' T energy-momentum flux through the 
closed hypersurface is 

6 



P T > 5 (E) +P S ' 1 (S 



s,r, 



1 1 dir ( — iA 2 Z)2£6T 

6 



Si 



<{-2pZ-^iZ+zA 2 Z)2(,Sy 



(5.28) 



. (5.29) 



5.3 Surface integral P 



5,6 



Thus we begin by calculating the integral P 5 ' 5 over the 3-tube S J , i.e., 



57T 



dir 



duj £ 2 F 5+ N T F 5 



(5.30) 



With the 5 parts of the field taken from d4.10M4TTTT) . and the normal to the 3-tube 
given by (14.121) . this becomes after some simplifications 

-,2 



P s ' s (E T ) = ^£dirljdu e ^(0 

x B(l -iu-Af{V- 2£aAa)B- ] 
which after some algebra gives 

xb(u- lioAo - 2£(z7- A)u + 4£ 2 (z7- A)aAa 



+ i 2 {V ■ AfV -2^{v ■ AfaAa)B 



6 



The angular integrals are then made using the formulas in Tabled] to get 



P 



n 



xB{ 0-2^A - 2^A + + iif^A 2 



2£ 3 —A 2 A)tp- 



(5.31) 



(5.32) 



(5.33) 
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which using (|4.1I) and (14.31) to revert from spinor into vector notation finally gives 



1 



1 



<l<T8\t)[--Z + ~iA 2 Z-±A 2 Z 



(5.34) 



We now repeat the same procedure for the integrals over the two 3-spheres S s 
closing the 3-tube at both ends, i.e., 

(5.35) 



With the 5 parts of the field taken from (I4.10M4~TTI) . and the normal to the 3-sphere 
given by (14.131) . this becomes after some simplifications 

r6 rr -2 



x B(l -£v-A) 2 {l -2i£aAa)B^ 



(5.36) 



which after some algebra gives 



P S > 5 {Z S ) 



xB 1 - 2i£aAa - 2£(v ■ A) + M^(u ■ A)aAa 



+ i\U-AY-2ii i (V-AYaAa)B 



6 



(5.37) 



The angular integrals are then made using the formulas in Tabled] to get 



dir - 2 S\0 



xB[ l-2i£~A + + Ai^iA 2 + t 2 ^A 2 - 2i^^A 2 A j B~ 



(5.38) 



which using (|4.1I) and (14.31) to revert from spinor into vector notation finally gives 



/ ^ <$2(f)( - ^riZ - - ^iA 2 Z 



1 

2^2- 



1 



i 

6' 



1 
10" 



T"2 



(5.39) 



Therefore, adding (15.341) and (15.391) . the P S ' S energy-momentum flux through 
the closed hypersurface is 

6 



P*' S (E) 



T2 



e" l dir , 

■ % 



1 Z + -iA 2 Z--^A 2 Z)eS 2 
3 10 ; ' 



6 



1 

6?" 



1 

6" 



z - —;z - -a 2 z - L iA ^z)e^ 



(5.40) 
(5.41) 
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5.4 Total surface integral 



Collecting the results from the previous subsections, i.e., Eqs. (15.151457T61) . (15.281 - 
I5T291) . and (l5T40ti54TT) . the three contributions P T ' T , P r ' 5 + P 5 ' r , and P 5 ' 5 to the 
surface integral P(X), i.e., Eq. (11.11) or the proper-time integral of the left-hand 
side of (11.51) . are: 



(5.42) 
(5.43) 



P T ' 5 (S) +P d ' 1 (S) =e 



<5,T/ 



2 dir {-Uj^zyi^i 



1 



1 



l 



(5.44) 
(5.45) 



P^E) = e 2 



dir 



1 
10 



6 



-iZ 



6 



10 ' 



T2 



(5.46) 
(5.47) 



When £ — ► the expression inside the parentheses in (15.421) is equal to (11.81) . i.e., 
the usual energy-momentum flow rate through a 3-tube surrounding the world-line 
flH[[3]|. This illustrates how much more complicated the rigorous formalism and 
calculations are than the customary ones. 



6 Volume integral: minus the r-integrated self-force 

In this section we calculate the two contributions P 5 ' r , and P s ' 5 to the volume 
integral (|1.4I) . i.e., minus the proper-time integral of the right-hand side of (11.51) . 

P(n) = dir JJJ^ d 3 tt^(JF - F + J). (6.1) 

However, since P = P + is bireal, we can anticipate that after integration the 
contributions arising from JF and — F + J will be equal. We therefore calculate 

P(Q) = J 2 dir J' d£ £ 2 U duo JF, (6.2) 
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where F is a sum of two terms as in (13.31) . and where the conserved current-density 
(13.21) is rewritten as a sum of three terms, i.e., 

J=4^{Ji + * + J*)> («-3) 
where, in quaternion 4-vector and spinor notations, using (14.11) . (14.31) . (14.91) . and 

(EH, 

j x = -(l - 2k£)Z£6' = -(1 - 2£z7 • X)BB + £5', (6.4) 
J 2 = i{£Z - Z\K)6 = B(i£A + 2ex°)B + 5, (6.5) 

J 3 = i£ K N T (^5 - -5') = A)B{V- 2£aAa)B + {5 - £5'). (6.6) 

Hence we calculate separate integrals of the form 

P« = e 2 r dir f'd£I%, (6.7) 
where II" is the angular integral 

K = ^-eII duJ ^ (6 " 8) 
in which F a is either F r or F 5 , i.e., 

F r = ^B*(v + 2£aAa)B + r, or F s = - £v-A)B*vB + 5. (6.9) 

Before starting to calculate the integrals (16.7H6.8I) let us illustrate our procedure 
by calculating the customary self-force in which J n reduces to the first term of 
(1641) . i.e., to J = -BB + £5', and the field F a to F r without the T factor, i.e., 
from ([O]) to F° = eB*(V + 2£aAa)B + /£ 2 . Referring to Tabled] for the angular 
integrals we get then 

n° = ^- JJ du {-BB + £5')^B*{v + 2£aAa)B + = ~\z5'. (6.10) 

Thus, interpreting <5'(£) as a Schwarz distribution and using (13.41) . we obtain 

/T2 l> + OC 1>T2 9 

rfir y d£U° = -j dire 2 -iZ\ (6.11) 

i.e., a self-force integral comprising the correct Schott term, but no Larmor term 
and no mass term as in (11.71) . Since the customary source current-density Jq = 
eZ5(^)/47r^ 2 corresponding to J by (16.31) does not contain the acceleration, it 
is clear that the self-force density JcF° cannot lead to a Larmor term, which is 
quadratic in the acceleration. This clearly demonstrates that the customary theory 
of point-charges in which the field is F° and the current density Jc is incomplete. 
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6.1 Angular integral 11^ 

From (f6~4j) and (l6\9l) . left, 

= JJ du{l- 2£v ■ A)BB + £5'^B*(v + 2£aAa)B + T 
= ~i ff duB(l-2£v-A)(v + 2£aAw)B + ^5'r 

= -B( + 2£-„4 - 2^-^* - oW Vt 
V 3 3 / £ 

= 0. (6.12) 



6.2 Angular integral II J 

From (1631) and (16~H . left, 

U " = hJJ duJ mA + 2 ^) B+ ^B*(u + 2^aAa)B + T 
= ^-JJdcu B(i£A + 2£ 2 xa)(v + 2£aAa)B + ^5r 
= i JJ duB ( l tAv + 2ii 2 AaAa + 2i?x^B + ^5i: 
= B(0 + 2teA^A)B + ^5T + ^- JJ dtu s(2^ 2 X a)i3+^T 
= -hA 2 ZST + JJ duo X K5T 

= -i(2A 2 Z + ^'ZjST (6.13) 
where we used 

j- JJ duo X K = -i(^A 2 Z + ^z). (6.14) 
To prove this formula, we start from (14.151) and (14.161) which immediately yield 

iff duJxK= i II duj { ~ A2 - ® ° ^ K 



-iA 2 Z - [[ du i{Z o S)S, (6.15) 



1 

4vr JJ 
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because integrating S = BuB + over solid angle gives zero. We therefore need to 
integrate 

i rr. ., _ iff 



— yy dw (z o s)u = — jj du §[zb*(-u)e\v 



JJ duJ (V[BZB*\ oV^V 



= \V[BZB*} = \{BZB* - B + ZB), (6.16) 
o o 

where we made a few elementary quaternion manipulations and, V[B'zB*} being 
independent on angular variables, used Tabled] to integrate. Thus 

i- JJdu(ZoS)S = ^(Z-ZZZ) = ^(Z + A 2 Z), (6.17) 

and combining with (16.151) we obtain (16.141) . 

6.3 Angular integral II J 

From (1631) and (|6"31) . left, 



nj = ^- ^ rfcj • .A)£(z7 - 2£aAa)B + {5 - £6')^B*{P + 2£aAa)B + T 



JJ duBi{v-A){v-2iaAa){v + 2iaAa)B + -^5-£5')T 
JJ duBi( K ~v-A-2ii{v-A)aAa-2ii{v-A)aAa^B + j{5-i5')T 
Bi(-0 - - 2i£-M 2 )s+~(<5 - £5')T 

~M 2 i(6 - £5') T - (6-18) 
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6.4 Angular integral 

From (ET4T) and (|6T9"1) . right, 



H s 1 = l- Jj du {1 - 2£P ■ A)BB + £5'hl - £P ■ A)B*PB + 6 
= Jl duB(l-2£v-A)(l-£v-A)vB + 5'5 
= JJ duoB(v-^(u- A)u + 2^(u- Afu^B+5'5 
= b(o-3£-A +0)b + 5'5 

= -CZ5'5 (6.19) 



6.5 Angular integral II^ 

From (1631) and d6H . right, 



= -— ^ duB(i£A + 2^ X a)B + 5-(l - ■ A)B*uB + 5 
= -±-JJduj B{iU + 2£ 2 X a)(l ~ iV-A)VB + Y 

= ~L II dU B {~ i ^ p + ^ 2 A{y ■ A)V - i 2 X 2ia + £ 3 (z7 ■ A)x2i^)B + ^5 2 
= B(-0+ i eA^A)B+±6 2 -^ JJ dujB(e{l-tf-A) X 2io)B + -^ 2 
= -i l -A 2 Ztb 2 - -L ff du (1 - iv- A) X Ki5 2 



= (iA 2 Z + -iZ - -iA 2 Z + • Z)Z^5 2 , (6.20) 
where we used (16.141) and the formula 



^ ft du{v-A) X K = - l -A 2 Z+ l -(z-Z)Z. (6.21) 
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To prove it, we start form (14.151) and (|4.16l) which yield 

1/7. Iff./ - . r \ 



— du X (v-A)K = — jj duo (-A 2 (u-A) -i'ZoS{V- A)^jK 

=iJJ^{-^-M 

+ ^ If duZo (S{u ■ Afj Z, (6.22) 



and directly leads to (16.211) because integrating [p ■ A)S = [p ■ A)BVB + over 
solid angle gives gives BAB + /3 = Z/3. 

6.6 Angular integral II3 

From (1631) and (I6T91) . right, 



U s 3 = -— jj dcu ■ A)B(v - 2£aAa)B + (5 - - A)B*uB + S 

= ~~L JJ dU Bl ( V -AW-2iaAa){l-iv-A)VB + i{5-£5')5 
= JJ duJ Bi{-V- A + £{u- A) 2 - 2i£(u- Apia 



+ 1ii 2 [V- A) 2 aAa)B + £{5 - £5')5 
-liA 2 Z + \tA 2 z)t{5-&')5. (6.23) 

6 5 / 



6.7 Total volume integral 

Combining the results of the above subsections, i.e., Eqs. (16.121) . (16.131) . (16.181) . 
(16.191 ). (16.201) . and (16.231) . the two contributions P T <5 , and P 5 ' 5 to the volume 
integral (11.41) . i.e., minus the proper-time integral of the right-hand side of (11.51) . 
are: 



PT2 1 C 

P T >\tt) = e 2 dir dn — iZ--iA 2 Z)5T (6.24) 

/T"2 /* ^2 1 

rfir y df(--U 2 z)f5'T (6.25) 
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and 



P 5 ' S (Q) = e 2 [ 2 dir [** dd--Z + -iA 2 Z - ^A 2 Z)£ 2 55' (6.26) 
J T1 v ? 3 5 / 

+ e 2 dir / 2 dti(+Uz + ^U 2 Z - ^A 2 Z + |(2 o I)z)^ 2 . 

(6.27) 



These integrals are rather complicated, what is of course necessary to match the 
complexity of the surface integrals. In fact, this could be expected from the 
conclusion of Ref. HU, where is was explained that the self-force density on an 
arbitrarily moving point-charge is given by \ (F + J — JF ) where F is the complete 
retarded field (12.71) . and where J is the full locally-conserved current density (12.81) 
or <E2k. 



7 The general energy-momentum flow 

In this section we transform the surface integrals (15 .42145471) and volume integrals 
(16.24j - [6~27l) according to the methods sketched in Sec.[3]and show that the sums 

P(E) = P T ' T (S) + P T ' 5 (£) + P 5 ' T (S) + P 5 ' 5 (S), (7.1) 
P(0) = P r '\tt) + P S >\Q), (7.2) 

are equal, i.e., P(X) = P(fi) := P as required by Ql.lb . and can both be put in 
the form: 



P(ti,t 2 ,6,6) 




(7.6) 

where all kinematical quantities are evaluated at the retarded proper time r r = r — £. 

For convenience we refer to (|7 .3147 .61) as the 'general energy-momentum flow' 
even though we could change its overall sign and call it the 'general r-integrated 



26 



self-force.' By 'general,' we mean (i) that this flow is a function of the parameters 
r i) r 2, £i> and £2! and (ii) that the function T, from which 5 derives by differentia- 
tion, belongs to a class of the nonlinear generalized function which may be further 
specified by physical conditions leading, for instance, to additional constraints on 
the Colombeau mollifier 77. 



7.1 Transforming the surface integrals 

We begin by transforming (15.421) into a double integral over £ and r by integrating 
by part. It easily comes 



dir 



2£' 



-iA Z)T 



T'2 



dir 



+ 



~2 



-2y 



dir 



6 



6 



dt[^Z--iA 2 Z){T 
di[^Z- 2 -zA 2 Z)r\ (7.7) 



We do the same with (15.431) and get, using (13.41 



«2 



6^ 



—Z + —iZ + -A 2 Z)T 

2£ 2 2£ 3 



+ r*jC«(-5(f ) ' + i ( ^) T 



n 



f2 



^ r y e _ 3 m2z j t2 - (v - 8) 

Therefore, adding (17.71) and (17.81) . the surface integrals (15.42145 431) can be written 

r & / 1 ~ 2 



P 



'/'> / <k(+j-Z--iA 2 Z)2<)T (7.9) 



We now rewrite (15.451) also as a double integral over £ and r so that (15.44145451) 
become 



P T ' 5 (E) = e 2 / dir / d£ 



T-2 



£2 



2e 



+ e z 



Tl 

ti V 3 



iZ + |U 2 ^)-)25T (7.10) 
/ 

(7.11) 
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Finally we rewrite (15.471) also as a double integral over £ and r so that (|5 .461 - 15^471) 
become 



ft 

(7.12) 

ft 



(7.13) 

Then, adding all contributions, i.e., (17.91) to (17.131) . we get P(£) in the form 

(EH. 

However, before verifying that the total volume integral can also be put in the 
form d7.3H7.6l) . an important remark is in order. Indeed, when adding (T7.9I) and 
(17. 101) a remarkable cancellation takes place: the l/2£ diverging terms cancel each 
other. This means that the Coulomb self-energy does not contribute to the total 
surface integral P(£). Indeed, the diverging term which will be interpreted as the 
'self-mass' will not be the usual Coulomb self-energy, i.e., the first term in (15.42b . 
but the leading divergence in (15.471) . which is an integral over a <5 2 -function. This is 
consistent with what was found in Ref . [6] for a point-charge at rest, namely that its 
self-energy is not the usual Coulomb self-energy but an integral over a 5 2 -function, 
so that the self-mass is fully concentrated at the position of the singularity. 



7.2 Transforming the volume integrals 

To show that d6.24H7q2"71) are equal to d7.3H7.6l) we need to transform them so that 
they do not contain 5' any more. We begin with (16.261) which, using (5 2 )' = 255' , 
can be integrated by parts as 



dir 



n 



6 



di [--Z+-iA z Z 



^A 2 z)e^' 



+ 



+ 



Tl 
T2 



Tl 

T> 



10 



<lir{—^Z+-iA 2 Z 

* 2 / 1 •■ 1 ... 



A'Z)Z 2 6 2 



dir 
dir 



ft 



10 



10 



{A 2 zy + ^-A 2 z - ^i{A 2 zy )i 2 5 2 



ft 



i • 



d£[+-Z--iA 2 Z + ^A 2 Z)£5 



(7.14) 
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Therefore, using (A 2 )- = 2Z oZ, Eqs. (l6T26l - t6?27l) become 



P s > s (£l) = e 2 



dir 



Tl 



1 



2C's iA2Z -To A2Z]e52 
1 ... 1 



dir 



(7.15) 



ft 



dn —z - -iz - -%A 2 z - j;(A 2 zy - ^-i(A 2 zy )£5 



2i 



6 



6 



10 



(7.16) 



For (16.251) a direct integration by part is not possible. However, if Eqs. (16 .24146351) 
and Eqs. (17. 151477161) are added together it is readily seen that their sum is equal to 
Eqs. (I7.3H7.6I) . This is because the identity 



(f <ST)' = 5T + £5'T + £8 2 



(7.17) 



and the integral 



ft 



ft 



%A 2 Z) (£<5T)' 



2 d$ (^ 2 z) (e*T) 



-\-iA 2 Z){^T) 



(7.18) 



enable to eliminate (|6.25l) . Thus, we have proved that P(O) = P T ' 5 (fi) + 
P 6 < 6 (£1) = P given by EM. 



8 The general self-interaction force 



By differentiating P{t\, r 2 , £i, £2) with respect to the proper time, that is by sup- 
pressing the ir integrations in (I7.3H7.6I) . we get the equality (11.51) . i.e., -P(S) = 
P(O) = -Q with 



Q(£i,6) = e 2 (ia 2 z)^r 



2 3 

ft ,1 



10 



4 



+ e 



ft 

2 d z (-^ z + + ^zy +y{A 2 zyy, 



(8.1) 
(8.2) 
(8.3) 
(8.4) 



where all kinematical quantities are evaluated at the retarded proper time r r = r — £. 
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Equation (I8.H - I8.4I) is the main result of this paper: it is an exact closed-form 
expression giving, in the limits £ 2 — > oo and £1 — > 0, the total self-interaction 
force on an arbitrarily moving relativistic point-charge. Because it leads to the 
Abraham-von Laue self-interaction force when taking these limits we call it the 
'general self-interaction force.' 

To calculate the limits £2 — » 00 and £1 — > in (l8.lH8.4t we first remark 
that if Z(r) is indefinitely differentiable in r, which implies that the velocity, the 
acceleration, and the biacceleration are smooth functions of r, the kinematical 
expressions inside the big parentheses are all smooth as functions of r r , that is as 
functions of £ as well. Then, because 



£5T 



0, and £5 



= 0, (8.5) 

the boundary terms d8.lH8.2l) give zero. As for (I8.3l - 48.4l) we use the formulas [UlEl 

f°° 1 

/ dr T(r)5{r)T(r) = -T(0), (8.6) 
Jo 2 

f°° dr 5 2 (r)T(r) = limC^] — + C[i]T'(0), (8.7) 
Jo e 

where C[o] and C[i] are functions of the Colombeau mollifier 77, i.e., 

f+00 r+00 
C[ ] = / dx7] 2 (—x), and C[i] = / dxxrf(—x). (8.8) 



In d8.6H8.7l) the function T(r) is a smooth function with compact support, which 
implies that T(±oo) = 0, what is the case of the kinematical expressions in 08.34 - 
8.41) if we suppose that Z(±oo) = 0, i.e., that the velocity is zero at r = ±00. 

Applying dOHO) to d8~TT - [8~4l) we finally get for the total self-interaction 
force Q self := Q(0, +00) the expression 

Q««(r) = e 2 (-lirn^Z(r) + ^iZ{r) + Uz{r) + ^A 2 (r)Z(r)), (8.9) 

where all kinematical quantities are evaluated at £ = 0, i.e., on the world-line 
at the proper- time r, and where we used (13.41) to replace the ^-differentiation in 
(18.71) by a ir-differentiation to get the term -\-\iC\^'Z{t). We therefore see that 
the <5 2 -integral (18.41) contributes two terms to the self-interaction-force (18.9b : a 
diverging contribution which can be interpreted as a self-mass as in the customary 
calculations, and a finite additional contribution to the Schott term which does 
not arise in the customary calculations. If these contributions do not need to be 
explicitly shown, one can rewrite (18.91) using (18.71) in the more concise form 



Q-*(r) = e 2 (-- ^ 6 2 (OZ(r r ) rf£ + -zZ(r) + -zA 2 (r)Z(r)). (8.10) 
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which like (I7.3H7.6I) is valid in any algebra of nonlinear generalized functions. 



9 The Lorentz-Dirac equation of motion 

An immediate application of the general formalism developed in this paper and 
of the self-interaction force Q self is the derivation of the equation of motion of a 
point-particle in arbitrary motion in an external electromagnetic field F ext . Indeed, 
the total force on a point-particle under the effect of its own field (12.71) and of an 
external field is 

m Z(r) = JJJ d 3 Q l - ((F + F ext ) + J - J(F + F ext )) 

= Q sdf (r) + ~e(F+(Z)Z(r) - Z{r)F^{Z)) , (9.1) 
with the subsidiary condition 

ZZ = 1, (9.2) 

which specifies that the world-line Z(r) is the world-line of a point in Minkowski 
space. In writing (19.11) we postulated that the total force could be written m§Z 
where m is a constant called the 'non-renormalized' mass, and that the external 
field F ext (X) is a smooth function of X — or at least smooth at any point X near or 
on the world-line Z{r) — so that the volume integral of the external force density, 
i.e., |(F+ t J - JF ext ) with J given by $2M or (Q, reduced to |e(F e + Z - ZF ext ) 
after an elementary integration over a three-dimensional 5-function. 

Substituting Q self from (18.91) we get therefore 

(m + Jim ^1) Z = e 2 (^kz + \%Z + \^ 2 z) + \e(F+Z - ZF ext ), (9.3) 

where we have put the diverging self-mass term on the left-hand side. We now 
use the subsidiary condition condition (19.21) . which by differentiation leads to the 
identities Z o Z = and Z o Z = —A 2 . Thus, multiplying (19.31 ) by Z and taking 
the scalar part we immediately obtain 

Sil + - = -. (9.4) 

2 6 3 

This equation is satisfied if the Colombeau mollifier i] is such that the moment 

r+oo 

Cm = / dxxri 2 (-x) = 1, (9.5) 
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which is an additional condition supplementing the minimal conditions (13.11) . 
Consequently, Eq. (19.31 ) becomes 

(m + lim -^)Z= \ie 2 (z + A 2 Z) + \e{F+ t Z - ZF ext ) , (9.6) 
v e^o ze / 6 z 

which after substituting the 'renormalized' or 'experimental' mass 

m := m n + lim — — , (9.7) 

is known as the Lorentz-Dirac equation of motion, i.e., Eq. (11.101) . where everything 
is evaluated at the proper time r. 

To fully appreciate what has been achieved, let us compare our derivation to a 
particularly elegant approach due to Barut |[T4|. The idea is to calculate the self- 
interaction force by analytically continuing the Lorentz force of a point-charge on 
itself, i.e., \e{F^ xf Z — ZFe\i) where F sel f is the customary Lienard-Wiechert field, 
to a point X = Z(t + e) close to the world-line. If this done, one easily finds [fT4l 
Eq.8] 

QtUr) = e 2 (-limiz(r) + Uz\r) + \iA\t)Z{t)), (9.8) 

which essentially differs from (|8.9I) by the absence of a finite contribution equiv- 
alent to |C[i]iZ(r) so that there is a problem in getting the correct Schott term. 
In fact, comparing with (18.101) . we see that the fundamental difference is that in 
Barut's approach the divering mass term is the 'Coulomb self-energy,' whereas in 
our rigorous formalism it is a '<5 2 -integral.' A similar situation arises in Dirac's 
approach [|4]|, in which instead of calculating the self-interaction force the flow 
of energy-momentum out of a 3-surface surrounding the world-line is calculated, 
i.e., Eq. (11.81) . which implies 

QllJj) = e 2 (- Jim iz(r) + 2 -zA\r)Z(r)) , (9.9) 

in which there is no Z{t) contribution at all. 

Therefore, with self-interaction forces like (19.81) or (19.91) . there is no other 
possibility than to make one or more additional suppositions in order to get the 
Lorentz-Dirac equation of motion. In a way or another this amounts to modifying 
or supplementing Maxwell's theory and/or the notion of causality which leads to 
the retarded potentials, what is not necessary with the self-interaction force (18.101) 
in which the 5 2 -integral supplies the missing Z(r) contribution. 
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10 Conclusion 



The results obtained in this paper show that the classical electrodynamics of point- 
charges is mathematically an internally-consistent theory, provided all potentials, 
fields, and currents are defined as nonlinear generalized functions. 

However, from a physical perspective, a mathematical point-charge, i.e., a 
simple pole like the potential (12.51) . is an idealization: a classical point-electron 
has an electric charge e, an intrinsic magnetic dipole moment ft, and spin 0. Its 
4-potential, expressed in the spinor notation of Sec.|4] is fT5l 

(10.1) 

T=T r 

It is therefore perhaps important to suggest that the methods used in this paper 
should be applied to the potential (110.11) . This would lead to an equation of 
motion generalizing that of Lorentz and Dirac, which could provide a physically 
satisfactory equation of motion for a classical electron with magnetic moment and 
spin. 
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Angular integrals 
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Table 1 : Angular integrals of expressions containing products of the idempotent 
a = \ ( 1 + iv) , the unit vector u(9, 0) , and of the vector A which does not depend 
on the angles 9 and 0. All expressions in which there is an odd number of z/factors 
are zero. Further identities may be obtained by taking the quaternion conjugate, 
complex conjugate, or quaternion biconjugate of X or Y and their integrals. 
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